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1. INTRODUCTION 


FEMOM3DR is a computer code written in FORTRAN 77 to compute 
electromagnetic(EM) radiation characteristics of antennas on a three dimensional object with 
complex materials (fig. 1) using combined Finite Element Method (FEM)ZMethod of 
Moments (MoM) technique[l]. This code uses the tetrahedral elements, with vector edge 
basis functions for FEM in the volume of the cavity and the triangular elements with the basis 
functions similar to that described in [2], for MoM at the outer boundary. By virtue of FEM, 
this code can handle any arbitrarily shaped three-dimensional bodies filled with 
inhomogeneous lossy materials. The basic theory implemented in the code is given in 
Appendix 1. 

The User’s Manual is written to make the user acquainted with the operation of the 
code. The user is assumed to be familiar with the FORTRAN 77 language and the operating 
environment of the computers on which the code is intended to run. The organization of the 
manual is as follows. Section 1 is the introduction. Section 2 explains the installation 
requirements. The operation of the code is given in detail in Section 3. Three example runs, 
the first EM radiation characteristics of an open coaxial line in a 3D PEC body, the second 
radiation characteristics of an open rectangular waveguide in a 3D PEC body, and the third 
EM radiation characteristics of an open circular waveguide in a three dimensional cavity are 
demonstrated in Section 4. Users are encouraged to try these cases to get themselves 
acquainted with the code. 


2. INSTALLATION OF THE CODE 

The distribution disk of FEMOM3DR is 3.5" floppy disk formatted for IBM 
compatible PCs. It contains a file named f emom3dr . tar . gz. This file has to be transferred 
to any UNIX machine via ftp using binary mode. On the UNIX machine, use the following 
commands to get all the files. 

gunzip f emom3dr . tar . gz 
tar -xvf femom3dr.tar 

This creates a directory FEM0M3DR-1 . 0, which in turn contains the 
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Figure 1 Illustration of cavity-backed radiating aperture in a 3D PEC body. 

The cavity is fed by a coaxial line or a rectangular waveguide or 
a circular waveguide at the input plane S inp . The fictitious outer surface 
S 0 is used to terminate the FEM computational domain. 
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subdirectories, FEM0M3DR (source files for the main code), PRE_FEMOM3 DR (source files 
for preprocessing code), Examplel, Example2 and Example3 As the code is written 
in FORTRAN 77, with no particular computer in mind, the source code in these directories 
should compile on any computer architecture without any problem. The code was successfully 
complied on a SGI machine, and the compilation can be done by using a makefile file for 
the different machines such as SUN, DEC or CONVEX etc. The complete listing of the 
directories in the distribution disk is given in Appendix 2. 

3. OPERATION OF THE CODE 

The computation of EM radiation characteristics from a specific geometry with 
FEMOM3DR is a multi-stage process as illustrated in figure 2. The geometry of the problem 
has to be constructed with the help of any commercial Computer Aided Design (CAD) 
package. In our case, we used COSMOS/M[3] as our geometry modeler and meshing tool. As 
FEMOM3DR uses edge based basis functions, the nodal information supplied by most of the 
meshing routines cannot be readily used. Hence, a preprocessor PRE_FEMOM3DR is written 
to convert the nodal based data into edge based data and then is given as input to 
FEMOM3DR. For the convenience of the users, who use different CAD/meshing packages 
other than COSMOS/M, PRE_FEMOM3DR accepts the nodal based data in a generic format 
also. The procedures involved for using COSMOS/M input data file or generic input data file 
are explained below. 

With the help of COSMOS/M, the geometry is constructed and meshed with 
tetrahedral elements. The user is assumed to be familiar with COSMOS/M package and its 
features. Once the mesh is generated, one needs to identify the following to impose proper 
boundary conditions: 

(a) tetrahedral elements with different material parameters 1 , 

(b) elements on PEC surfaces 

(c) elements on the outer boundary 

(d) elements on the input plane 

1. COSMOS/M has a feature by which it can group tetrahedral elements with different material proper- 
ties into different groups. For a generic file input, the user has to specify the material property index for 
each tetrahedral element to indicate its material property group(see Appendix 4). 
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Generic Mesh data ^ 

1 from any CAD/ 

I Meshing Program 
\ * J>IN file ^ 

Figure 2 Flow chart showing the various steps involved in using FEMOM3DR 

This is done using the available features in COSMOS/M. Sample *.SES files of 
COSMOS/M which illustrate these features are given in Appendix 3. Finally the *.MOD file 
is generated with the required mesh information. PRE__FEMOM3DR accepts the *.MOD file 
as input and generates the required edge based data. 

For users, who can do geometry modelling and meshing of the model with any other 
CAD package, the nodal based information is required to be placed in a file pnoblem. PIN, 
where problem is the name of the problem under consideration. The format required for 

. PIN file is given in Appendix 4. Note that all the dimensions of the geometry are assumed 
to be in centimeters. 

The PRE_FEMOM3DR code gives the following prompts: 

pre_f emom3 dr 

Give the problem name: 

The problem name is the user defined name for the particular problem under consideration. 

COSMOS file (1) or GENERIC (2) file? 

If you are using *.MOD file from COSMOS/M, give 1 or using the generic input data file 
explained above, give 2. 

TFMOM3 DR generates the following files with required edge based information. 
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(a) probl em_n.od.al . DAT - Node coordinates and the node numbers for each element 

(b) problem_edges . DAT - Information on edges, such as nodes connecting each edge, etc. 

(c) probl em_suz£ ed . DAT - Information on edges on outer surfaces. 

(d) probl em_surf el . DAT - Information on edges on input surface. 

(e) probl em.POUT - General information on the mesh. 

The files (a) to (d) are used as input for FEMOM3DR. Users need not interact or modify the 
above files. 

After PRE_FEMOM3DR is run, all but one input data file required for FEMOM3DR 
are ready. FEMOM3DR expects to find problem . MAT file which contains the material 
constants information required for the volume elements. The format of the probl emMAT is 
as given below: 

N g , Maximum number of material groups 

*vi > Ur i Complex relative permittivity, complex relative permeability respectively 

e r2’ Pr 2 for material groups 1, 2, 3, , N n 

8 

^ rN g ’ PrN g 


In the PRE_FEMOM3DR, all the terahedral elements are given the material group index. The 
material parameters given in problem. MAT are read into FEMOM3DR and the proper 
material parameters are assigned to each tetrahedral element according to its material property 
index. Once the problem. MAT is ready, FEMOM3DR code can be run. The FEMOM3DR 
code gives the following prompts: 

f emom3dr 

Give the problem name : 

This name should be the same as given for PRE_FEMOM3DR 

Frequency (GHZ) : 

This is the frequency of operation. If the dimensions of the problem are in wavelengths, 
frequency should be specified as 30 GHz as FEMOM3DR assumes that all dimensions are in 
centimeters. 
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Give the type of feed line : 
coax ( 1 ) , rect wg(2), cir wg(3) 

This is to specify the type of feed line to be used. User should give 1 if coaxial feed is used, or 
2 if rectangular waveguide is used as feed, or 3 if circular waveguide is used as feed. Depending 
on the feed line to be used, FEMOM3DR gives different prompts to input the feed line 
parameters. 

For coax(l) 


Coaxial feed line 

Give Inner rad, rl(cm). Outer rad, r2 (cm) : 

Specify the inner radius and outer radius of the coaxial line. 
Dielectric const for the coaxial line, erl 
Specify the dielectric constant used for the coaxial line. 


For rect wg(2) 

Rectangular waveguide feed 

Give waveguide dimensions : a(cm), b(cm) 

Specify the waveguide dimensions, broad wall dimension a(cm), narrow wall 
dimension b(cm) 

For cir wg(3) 

Circular waveguide feed 

Give the radius of circular waveguide aa(cm): 

Specify the radius of the circular waveguide in cms. 

For computing radiation pattern, give Theta (degs) - 
start angle, stop angle and increment 

Specify the start and stop angles of 0 in degrees. Radiation patterns will be computed in both 
<!> = 0 and <j) = 90° planes. 

FEMOM3DR generates the file problem. OUT, which contains information on 
CPU times for matrix generation, matrix fill, the input characteristics and the radiation pattern 
data. FEMOM3DR also generates another file probIem_bicgd.DAT which contains 
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information on convergence history of diagonally preconditioned biconjugate gradient 
algorithm used to solve the matrix equations. 

4. SAMPLE RUNS 


Three example runs are illustrated in this section. They are selected to illustrate 
some of the features of FEMOM3DR. 

Example 1 : Radiation from an open coaxial line in a 3D PEC body 


z 



OUTER BOUNDARY FOR FEM-MoM : 4.5cmX4.5cmX 1.5cm 


Figure 3 Open coaxial line in a 3D PEC box. The analysis is carrired out at 5.73GHz. 
Inner radius of the coaxial line is 1.0cm and outer radius is 1.57cm. 


An open coaxial line in a 3D PEC box is considered. Assuming the dominant TEM 
mode propagation in the coaxial line the radiation pattern and input characteristics (z=0 as 
reference plane) are calculated. 
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First the PRE_FEMOM3DR 

c j r@caph : { 53 } pre_f emom3dr 
Give the problem name : 
coax 

COSMOS file(l) or GENERIC (2) file ? : 

1 

Opening file : coax. MOD 

Read the following data 
Nodes= 403 

Elements= 1201 

Elements on surface 1 ~ 468 

Elements on surface 2= 44 

Max number of material groups^ 1 

Forming the edges ! ! ! Be patient ! ] ! 

Order of the FEM matrix- nptrx= 1554 

************************************************** 


Number of nodes = 403 

Number of elements^ 1201 

Number of total edges= 2001 

Number of elements on Surface 1= 468 

Number of elements on Surface 2= 44 

Number of edges on surface 0(pec)= 447 

Number of edges on surface 1= 702 

Number of edges on surface 2 - 84 

Max number of material groups = 1 

Order of FEM matrix^ 1554 

Order of MoM matrix (electric cuurent ) = 702 

Unknown for the magnetic current^ 702 

Number of unknowns on Input plane= 48 

Order of Hybrid FEM/MoM matrix^ 2256 

************************************************** 


The coax . MAT file for this problem is given below: 

1 

( 1 . 0 , 0 . 0 ) ( 1 . 0 , 0 . 0 ) 
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And then FEMOM3DR : 

cjrOcaph: {55} femom3dr 

Give the problem name : 
coax 

Reading the input ! ! 

Finished reading the data 
Give frequency of operation : GHz 
5.73 

Give the type of feed line : 
coax(l) , rect wg(2), cir wg(3) 

1 

Coaxial feed line 

Give Inner rad, rl(cm). Outer rad, r2 (cm) : 

1.0 1.57 

Dielectric const for the coaxial line, erl 

1.0 

For Computing the radiation pattern, give Theta(degs)- 
Give start angle, stop angle and increment 

-180 180 10 

*************************************** 

* 

* 

* 

: k 

* FEMoM3DR (Version 1.0) * 

Problem : coax 

* ( BiCGDNS Solver) * 

* * 

* * 

*************************************** 

RADIATION CHARACTERISTICS OF AN ANTENNA ON 
A 3D BODY USING FEM/MOM HYBRID METHOD 

Frequency (GHz) = 5.730000 

Order of the FEM-MoM matrix^ 2256 

Order of the MoM matrix = 702 

Coax feed is used 

with characteristic impedance (ohms ) = 27.06454 

Radius of inner conductor (cm) = 1.000000 

Radius of outer conductor (cm) = 1.570000 

Dielectric constant = 1.000000 
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Generating FEM matrix 

Number of non zeros in amat ( zmatrices ) = 

Time to fill FEM matrix ( secs ) = 0.2383671 

net= 1554 

Time to fill zmatrixeh= 5 . 1747322E-02 
Generating Zmatrices 
Entering zmatrixej 

Time to fill zmatrixej ( secs ) = 155.7325 

Entering zmatrixem 

Time to fill zmatrixem ( secs ) = 221.0665 

Time to fill zmatrices ( secs )= 377.7292 

Total no of non zeros after adding zmatrices^ 900366 

Calling selmts_coax 
Entered selmts.f 

betal0= 1.200088 r2= 1.570000 

rl= 1.000000 zc= 27.06454 

Out of selmts 

Total nonzeros in amat after smat= 902502 

Solving the system of equations Ax=B by BiCGDNS 
CONVERGENCE ACHIEVED in 1833 iterations 

Residual Norm= 6 . 8373792E-04 
Time to solve by BiCGDNS ( secs ) = 1133.188 

Input parameters for the coaxial feed 

Reflection Coefficient Sll= (0.2084835,-0.6935343) 

Return Loss (db) = -2.802917 

Normalized Input Admittance , Yin/Yo= (0.2449467,0.7144601) 
Normalized Input Impedance, Zin/Zo= (0.4293905,-1.252445) 


RADIATION PATTERN (phi=0 deg plane) 
Theta (deg) 101og|Eth| ^2 101og|Eph| A 2 


-180 

-54.19743 

-48.60585 

-170 

-10.79702 

-46.29274 

-160 

-5.749784 

-45.47018 

-150 

-3.863171 

-46.03180 

-140 

-3.623751 

-48 . 02142 

-130 

-4.477120 

-51.30671 

-120 

-5.877765 

-52.62412 

-110 

-6.739963 

-49 . 10905 

-100 

-5.970154 

-45.52240 

-90 

-4.136563 

-42 . 71252 
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-80 

-2.249640 

-40.37156 

-70 

-0.6981635 

-38.32143 

-60 

0.4049507 

-36.53478 

-50 

0 . 9726373 

-35.07118 

-40 

0.8663030 

-34 . 02491 

-30 

-0.1586271 

-33.50406 

-20 

-2 . 607281 

-33 . 62665 

-10 

-8 . 010056 

-34.52869 

0 

-36.21795 

-36.36510 

10 

-7 . 811026 

-39.21128 

20 

-2 . 518029 

-42.22832 

30 

-0.1083566 

-42.46354 

40 

0 .8974382 

-40 . 67097 

50 

0.9943753 

-39.39063 

60 

0.4234217 

-38 . 98084 

70 

-0 . 6786989 

-39.25732 

80 

-2 .226301 

-40.00126 

90 

-4.108023 

-41 .00850 

100 

-5.937835 

-42.10268 

110 

-6 .708505 

-43.18414 

120 

-5 . 847709 

-44.28653 

130 

-4 . 445064 

-45.59224 

140 

-3 .588042 

-47.41969 

150 

-3 . 822448 

-50 . 17274 

160 

-5 . 698858 

-53.34728 

170 

-10.70996 

-52.22258 

180 

-54.19810 

-48.60595 


RADIATION PATTERN 

(phi=90 deg pi 

(deg) 

lOlog | Eth | A 2 

10 log | Eph | A 

-180 

-48 . 60547 

-54.19707 

-170 

-10 . 82473 

-48.70198 

-160 

-5.750631 

-46.34335 

-150 

-3.852743 

-46.03323 

-140 

-3 . 605450 

-47.45369 

-130 

-4 . 451667 

-49.95887 

-120 

-5 . 845438 

-50.09461 

-110 

-6.704795 

-46.91512 

-100 

-5.939765 

-43.87301 

-90 

-4.112114 

-41.47113 

-80 

-2.229023 

-39.44118 
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lete session of this run on a SGI 




Example 2 : Radiation from an open rectangular waveguide in a 3D PEC 
body 


z 



0.9X 

OUTER BOUNDARY FOR FEM-MoM : 1.0AX0.6AX0.5X 


Figure 4 Open Rectangular waveguide in a 3D PEC box 


An open rectangular waveguide in a 3D PEC box is considered. Assuming the 
dominant TE 10 mode propagation in the waveguide the radiation pattern and input 
characteristics (z=0 as reference plane) are calculated. 


First the PRE_FEMOM3DR 

c j r@caph : { 11 } pre_f emom3dr 
Give the problem name : 
rwg 

COSMOS file(l) or GENERIC (2) file ? : 

1 

Opening file : rwg. MOD 

Read the following data 
Nodes= 565 
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Element s= 
Elements on 
Elements on 
Max number 


1743 

surface 1= 560 
surface 2= 56 
of material groups= 


Forming the edges ! ! ! Be patient ! ! ! 


1 


Order of the FEM matrix- nptrx= 2159 

************************************************** 


Number of nodes= 565 

Number of elements^ 1743 

Number of total edges= 2836 

Number of elements on Surface 1= 560 

Number of elements on Surface 2 = 56 

Number of edges on surface 0(pec)= 677 

Number of edges on surface 1= 840 

Number of edges on surface 2= 95 

Max number of material groups= 1 

Order of FEM matrix= 2159 

Order of MoM matrix (electric cuurent)= 840 

Unknown for the magnetic current = 840 

Number of unknowns on Input plane= 73 

Order of Hybrid FEM/MoM matrix^ 2999 

************************************************** 
cjr@caph: {12 } 


The rwg . MAT file for this problem is given below: 

1 

( 1 . 0 , 0 . 0 ) ( 1 . 0 , 0 . 0 ) 

And then FEMOM3DR : 

cjrdcaph: {18} femom3dr 
Give the problem name : 
rwg 

Reading the input ! ! 

Finished reading the data 
Give frequency of operation : GHz 
30.0 

Give the type of feed line : 
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coax ( 1 ) , rect wg(2), cir wg(3) 

2 

Rectangular waveguide feed 
Give waveguide dimensions : a(cm), b(cm) 

0.7 0.31 

For Computing the radiation pattern, give Theta (degs)- 
Give start angle, stop angle and increment 

-180 180 10 




FEMoM3 DR (Version 1.0) * 

Problem : rwg 

* ( BiCGDNS Solver) * 

* * 

* * 

*************************************** 

RADIATION CHARACTERISTICS OF AN ANTENNA ON 
A 3D BODY USING FEM/MOM HYBRID METHOD 

Frequency (GHz) = 30.00000 

Order of the FEM-MoM matrix^ 2999 

Order of the MoM matrix = 840 

Rect w/g feed is used 

a(cm)= 0.7000000 b(cm)= 0.3100000 


* 


* 


Generating FEM matrix 

Number of non zeros in amat ( zmatrices ) = 
Time to fill FEM matrix ( secs ) = 0.3556371 

net= 2159 


Time to fill zmatrixeh ( secs ) = 
Generating Zmatrices 
Entering zmatrixej 
Time to fill zmatrixej ( secs ) = 
Entering zmatrixem 
Time to fill zmatrixem ( secs ) = 
Time to fill zmatrices (secs) = 


6 . 23 65055E-02 

219.9152 

316.8897 

537.7435 


Total no of non zeros after adding zmatrices= 
calling selmts_rwg 
Entered selmts.f 

Total nonzeros in amat after smat= 1398430 


25494 


1393424 
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Solving the system of equations Ax=B by BiCGDNS 
CONVERGENCE ACHIEVED in 1172 iterations 

Residual Norm= 8 . 4834168E-04 
Time to solve by BiCGDNS= 1134.912 

Input parameters for the Rect W/G feed 

Reflection Coefficient Sll= (-9 . 9301338E-05, -0 .2776211) 
Return Loss (db) = -11.13095 

Normalized Input Admittance , Yin/Yo= (0.8570415,0.5156051) 
Normalized Input Impedance , Zin/Zo= (0.8567256,-0.5154150) 


RADIATION PATTERN (phi=0 deg plane) 


Theta (deg) 

-180 

-170 

-160 

-150 

-140 

-130 

-120 

-110 

-100 

-90 

-80 

-70 

-60 

-50 

-40 

-30 

-20 

-10 

0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 


lOlog | Eth | A 2 

-61.92020 

-61.31005 

-63 .44328 

-68.74234 

-70.03734 

-66.70441 

-66.59460 

-69.86186 

-80.20856 

-75 . 01617 

-66.67963 

-61.93605 

-58 . 61893 

-56.34359 

-54.92506 

-54.19392 

-54 . 09115 

-54.73644 

-56.31515 

-58 . 89853 

-62.37308 

-66.48820 

-68.71347 

-65.26954 

-61 . 62342 

-59.33685 

-58.21431 

-58.07228 

-58 . 69095 


lOlog | Eph | A 2 

-15.32931 
-15 . 91580 
-17 . 62041 
-20.02731 
-21.95038 
-22.52033 
-22.63553 
-22.92205 
-22.88857 
-21 . 71360 
-19.52076 
-17 . 02657 
-14 . 62941 
-12.44571 
-10.51392 
-8.882670 
-7 . 628004 
-6.834783 
-6 . 567 93 8 
-6.851701 
-7 . 662488 
-8.935760 
-10.58671 
-12 . 53923 
-14 . 74527 
-17 . 16811 
-19.69254 
-21 . 91248 
-23 . 08877 
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-19.29325 

-63.71844 

140 

-23 . 87317 

-61 . 69011 
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-22.24701 

-60.44720 

160 

-18.23417 

-60 . 01572 

170 

-16.01494 

-60.42999 

180 

-15.32931 

-61 . 92015 


* 


* 




The complete session of this run on a SGI machine along with all the files is kept 
in the directory . / FEMOM3 DR- 1 . 0/Example2 . 

Example 3: Radiation from an open circular waveguide in a 3D PEC box 


z 



OUTER BOUNDARY FOR FEM-MoM : 10cmX10cmX4cm 


Figure 5 An open circular waveguide in a 3D PEC box 

An open circular waveguide in a 3D PEC box is considered. Assuming the 
dominant TE n mode propagation in the waveguide the radiation pattern and input 
characteristics (z=0 as reference plane) are calculated at 2.8GHz. 
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First the PRE_FEMOM3DR 

c j r@caph : { 60 } pre_femom3dr 
Give the problem name : 
cwg 

COSMOS file(l) or GENERIC (2) file ? : 

1 

Opening file : cwg. MOD 

Read the following data 
Nodes= 601 

Elements^ 1915 

Elements on surface 1= 564 

Elements on surface 2 - 64 

Max number of material groups= 1 

Forming the edges ! ! ! Be patient ! ! ! 

Order of the FEM matrix- nptrx= 2332 

************************************************** 


Number 

of 

nodes= 

601 


Number 

of 

elements: 

= 1915 


Number 

of 

total edges= 3071 


Number 

of 

elements 

on Surface 1= 

564 

Number 

of 

elements 

on Surface 2 - 

64 

Number 

of 

edges on 

surf ace 0 (pec ) = 

739 

Number 

of 

edges on 

surface 1= 

846 

Number 

of 

edges on 

surface 2= 

106 


Max number of material groups = 1 

Order of FEM matrix^ 2332 

Order of MoM matrix ( electric cuurent)= 846 

Unknown for the magnetic current = 846 

Number of unknowns on Input plane= 86 

Order of Hybrid FEM/MoM matrix^ 3178 

************************************************** 
The cwg . MAT file for this problem is given below: 

1 

( 1 . 0 , 0 . 0 ) ( 1 . 0 , 0 . 0 ) 
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cjrOcaph: {65} femom3dr 
Give the problem name : 
cwg 

Reading the input ! J 
Finished reading the data 
Give frequency of operation : GHz 
2.8 

Give the type of feed line : 
coax ( 1 ) , rect wg(2), cir wg(3) 

3 

Circular waveguide feed 
Give the radius of circular waveguide aa(cm) : 

3.75 

For Computing the radiation pattern, give Theta (degs ) - 
Give start angle, stop angle and increment 

-180 180 10 

*************************************** 


FEMoM3 DR (Version 1.0) * 

Problem : cwg 

* ( BiCGDNS Solver) * 

* * 

* * 

*************************************** 

RADIATION CHARACTERISTICS OF AN ANTENNA ON 
A 3D BODY USING FEM/MOM HYBRID METHOD 

Frequency (GHz) = 2.800000 

Order of the FEM-MoM matrix^ 3178 

Order of the MoM matrix = 846 

Circular w/g feed is used 
Radius of the w/g(cm)= 3.750000 


* 


* 


Generating FEM matrix 

Number of non zeros in amat ( zmatrices ) = 
Time to fill FEM matrix ( secs ) = 0.3919129 

net= 2332 


Time to fill zmatrixeh= 6 . 6025734E-02 
Generating Zmatrices 
Entering zmatrixej 


27988 
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CONVERGENCE ACHIEVED in 2245 iterations 

Residual Norm= 8 . 8973634E-04 

Time to solve by BiCGDNS ( secs ) = 2174.860 

Input parameters for the cir waveguide feed 

Reflection Coefficient Sll= ( -0 . 1300059 , 3 . 4798384E-02 ) 

Return Loss (db) = -17.42023 

Normalized Input Admittance, Yin/Yo= ( 1 . 295194 , -9 . 1804117E-02 ) 
Normalized Input Impedance , Zin/Zo= ( 0 . 7682255 , 5 . 4452278E-02 ) 


RADIATION PATTERN (phi=0 deg plane) 
Theta (deg) 101og|Eth| A 2 101og|Eph| A 2 


-180 

-22.75692 

-63.35085 

-170 

-24.54149 

-60.63559 

-160 

-30.61100 

-59.22809 

-150 

-27.98133 

-58.73494 

-140 

-21 . 97158 

-58 . 88988 

-130 

-19 . 15922 

-59.52450 

-120 

-17.95044 

-60.42453 

-110 

-17 . 54107 

-61 . 02032 

-100 

-17.34824 

-60.33348 

-90 

-16.88481 

-58.21341 

-80 

-15.86516 

-55 . 56404 

-70 

-14.29549 

-53.00042 

-60 

-12.37232 

-50.69931 

-50 

-10.33036 

-48 . 69147 

-40 

-8.376839 

-46 . 99461 

-30 

-6.681361 

-45.64932 

-20 

-5.375779 

-44.71483 

-10 

-4.553909 

-44.25065 

0 

-4.271783 

-44.29565 

10 

-4.548341 

-44.85369 

20 

-5.365640 

-45.88722 

30 

-6.668044 

-47.32177 

40 

-8.360756 

-49 . 06322 

50 

-10.30902 

-51.02592 

60 

-12.33875 

-53.16288 

70 

-14.23919 

-55.49183 

80 

-15.77886 

-58.11651 

90 

-16.77509 

-61.16561 

100 

-17.23692 

-64.18618 

110 

-17 . 45133 

-64 . 79224 
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120 

-17 . 89702 

-6 2.96498 

130 

-19.15180 

-61.55912 

140 

-22.02265 

-61.53593 

150 

-28.06933 

-63.19190 

160 

-30.24335 

-66.42630 

170 

-24.41655 

-67 . 00469 

180 

-22.75692 

-63.35098 


* 


* 


RADIATION PATTERN (phi=90 deg plane) 


Theta (deg) 

-180 

-170 

-160 

-150 

-140 

-130 

-120 

-110 

-100 

-90 

-80 

-70 

-60 

-50 

-40 

-30 

-20 

-10 

0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 


lOloglEthj ^2 

-63.35092 

-61.79300 

-62.38220 

-65.55884 

-64.24979 

-58.26767 

-54.45654 

-52.24393 

-51.02455 

-50.39310 

-50 .00290 

-49.51344 

-48.66899 

-47.46519 

-46.14532 

-45.00156 

-44.24279 

-43.98824 

-44.29570 

-45.18002 

-46.61389 

-48 . 51177 

-50 . 69865 

-52.88061 

-54 . 68048 

-55.80803 

-56.25292 

-56.26787 

-56.18223 

-56.26547 

-56.69032 


10 log | Eph | ^2 

-22.75693 

-23.50035 

-25.78485 

-29.34372 

-32 . 63723 

-33.54564 

-31.86797 

-27.87085 

-23 . 65279 

-20.07571 

-17.06390 

-14.43425 

-12.06949 

-9.930590 

-8 . 040731 

-6.461413 

-5.266761 

-4.522365 

-4.271783 

-4.530254 

-5.283057 

-6.487364 

-8 . 078716 

-9.984387 

-12.14472 

-14.53915 

-17.21146 

-20.28884 

-23.97681 

-28.39919 

-32.71435 
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130 

-57.53803 

-34.63153 

140 

-58.86083 

-33.68311 

150 

-60.81037 

-29.88833 

160 

-63.46891 

-26.00092 

170 

-6 5.06511 

-23.57710 

180 

-63.35083 

-22.75694 


* 


* 




The complete session of this run on a SGI machine along with all the files is kept in the 
directory . / FEMOM3 DR - 1 . 0/Example3 . 


5. CONCLUDING REMARKS 

The usage of FEMOM3DR code is demonstrated so that the user can get acquainted 
with the details of using the code with minimum possible effort. As no software can be bug 
free, FEMOM3DR is expected to have hidden bugs which can only be detected by the 
repeated use of the code for a variety of geometries. Any comments or bug reports should be 
sent to the author. As the reported bugs are fixed and more features added to the code, future 
versions will be released. Information on future versions of the code can be obtained from 

Electromagnetics Research Branch (MS 490) 

Flight Electronics and Technology Division 
NASA-Langley Research Center 
HAMPTON VA 23681 
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Appendix 1 


Theory for FEMOM3DR 

This appendix is intended to give a brief description of the theory behind the code. 
The geometry of the structure to be analyzed is shown in figure 1. S B represents the area of 
the fictitious outer boundary to be used for terminating the FEM computational domain and 
^ inp represents the area of the input plane. The electric field inside the compuational domain 
satisfies the vector wave equation [4] 

Vx (ir r VxE )-^ E = 0 a) 

where e r and are the relative permittivity and relative permeability of the medium. The 
time dependency of exp (j m) is assumed through out this report. To facilitate the suitable 
solution of the partial differential equation in (1) via FEM, multiply equation (1) with a vector 
testing function T and integrate over the volume of the computational domain. By applying 
suitable vector identities, equation(l) can be written in its weak form as, 

JJJ * (jf VxE^eJJ/T • Edv = x A Vxe)/v (2) 

Applying the divergence theorem to the right hand side of equation(2), the volume integral is 
written as sum of the surface integral over the surface S 0 terminating the FEM computational 
domain and the surface integral over S inp at the input plane. 



where h Q is the unit outward normal to the surface S 0 and h i is the unit outward normal to the 
surface S inp . 

To discretize the above volume and surface integrals, the FEM computational domain 
is subdivided into small volume tetrahedral elements. The electric field is expressed in terms 
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of vector edge basis functions [2] which enforce the divergenceless condition of the electric 
field implicitly 

6 

E = I «/W, (4) 

i = 1 

where e t ’s are the unknown coefficients associated with each edge of the tetrahedral element 
and W- ’s are the basis functions and are given in detail in [5]. The testing function T is taken 
to be the same set of basis functions as given in equation (4), i.e., 

T = W- j=l,2,3,4,5,6 (5) 

The discretization of the FEM computational volume automatically results in 
discretization of surfaces S 0 and S- np in triangular elements. The evaluation of the surface 
integral over the outer boundary is carried out either by using Method of Moments(MoM) and 
the evaluation of the surface integral over the input plane is carried out using mode matching 
method. 


Evaluation of surface integral over S Q - MoM formulation: 

At the fictitious outer boundary the electric field is subjected to the condition that the 
fields are continuous across the boundary, i.e., 

El = El (f,) 

I at S + a \at S' 0 ^ 

where S + Q denotes the outer side of S Q and S 0 denotes the inner side of S a . The electric field 
El is the field quantity being evaluated in the computational volume through FEM. The 
electric field ouside S Q is evaluated explicitly using the following equation[4, eq.3-83]: 


E| = - VxF - joou^A + — - — VV • A 

I at S ° /Vrm 




(7) 


where 


and 


A = Magnetic Vector Potential = —- J J 


1 f f Jexp(-yfc„|r-r 0 |) 


r- r. 


ds 


( 8 ) 
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F = Electric Vector Potential 




Ms\p(~ jk a 

1 

O 


r ~ r o 



ds 


( 9 ) 


J and M are assumed to be equivalent electric and magnetic currents respectively at the outer 
surface S Q . The equivalent currents radiating in free space are used in the equation (7) to 
compute the electric field outside V (figure 6). 

Substituting equation (7) into equation (6) and multiplying by a testing function 
n 0 x T on both sides and integrate over the surface S 0 , results in: 


J J (n„ x T ) • = -JJ(« 0 xT).(VxF)rfj-yo)(i J f(» xT)»A<& 

s. I 

+ ?irlJ ( ” oXT) * (VV * A) ^ < 10 ) 

J o So 

After some mathematical manipulations [6, pp.42], [7, pp.135], and substituting equations (8) 
and (9) in the above equation, it can be rewritten as: 



x 





{V * ( ^ xT)} {j/ 


ho 


(11) 


where 


ff 


indicates that the singular point has been removed and 


G = 


e*p(-# 0 |r-r 0 [) 


( 12 ) 


Equation (1 1) is written in a matrix form by choosing the proper basis functions for M and J 
and accordingly using the testing function h Q x T . Within each surface triangle, the surface 
currents can be expressed as 
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Fictitious outer boundary S Q 


Figure 6 Equivalent current representation of the outer surface S o 
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( 13 ) 


3 

M = Exh 0 = -£ e .'(«o xW i) 

i= 1 
3 

J = X hiK X w,.) ( 14 ) 

i = 1 

and the testing function as 

h 0 xT = n o x\V j j=l,2,3 (15) 

In equation (13), g. represents the same unknown coefficient as in equation (4) and in 
equation(14) /• represents the unknown coefficient for the surface electric current densisty. In 
equations (13) and (14), it is interesting to note that, the vector edge basis functions W-, 
which are initially used for electric field are used to represent the surface current densities in 
the form of h Q x W f . The expansion functions W . are used to build tangential continuity 
into the field representation. In contrast, the cross product of h 0 with these functions results 
in another set of basis functions which guarantee normal continuity with zero curl and 
nonzero divergence and hence are ideally suited for representing surface current densities [2]. 
During the current investigation, it has been observed that the roof top basis functions for 
triangular patches used by Rao[6] and the basis functions used here proved to be numerically 
identical to each other confirming the above point of view. 

Equations (13-14) are substituted in equation (1 1) and integrated over all the triangular 
patch elements on surface S Q to obtain the following matrix equation: 


where 


Mi 


{e} + 


M~ 


{/} = { 0 } 


(16) 



C ^ o s 


+ { v * ( "o x T »{/J (V-J)G* 1 j* 


(17) 


(18) 
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and {0} is the null vector. The singularities in evaluating the integrals in equation (18) are 
handled analytically by using the closed form expressions given in [8], 

Using Maxwell’s equation VxE = — 7'coq o jLl r H, the surface integral on the right hand 
side of the equation (3) can be written as 

“U T *(« 0 xivxE^'-J J T.{a,.xlv X E)* = J J T • (n o x H)efc 

Sinp Y S a 

+ J J T . x H inp )ds (19) 

Sinp 

where H inp is the magnetic field over the input plane obtained from matching the modal 
expansion of waveguide fields with the unknowns fields at the input plane[9]. By equivalence 
principle, it can be noted that J = h Q x H on the surface S a . Substituting this into equation 
(19), equation (3) can be rewritten as: 


JJj jJ-(VxT)*(VxE)</v-^e r JJjT»Edv = JJt.J^ + J J T • («,. x H fBp )* (20) 
v r V s„ s lnp 

Substituting equations (4), (5) and H inp in the above equation and integrating over all the 

tetrahedral elements to evaluate the volume integrals on the left hand side and integrating over 

all the surface triangular elements to evaluate the surface integrals on the right hand side, it 

can be written in a matrix form as 


F >] { 


e} + 


{/} = (M 


( 21 ) 


where [F j] includes the volume integration and the surface integration over the input place 
due to mode matching, 


[f 2] = Jl ' T * ids (22) 

and {&;} is the excitation vector due to the dominant mode incident in the waveguide. The 
evaluation of the volume integrals over a tetrahedral element is given in detail in [5]. 
Equations (21) and (16) are combined to form a system matrix equation: 


Mi 


M r 


0 

h 


( 23 ) 
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In the above system matrix F j and F 2 are sparse matrices and M j and M 2 are dense matrices 
and also the total matrix is complex and non- symmetric in nature. This matrix equation is 
solved using a diagonally preconditioned biconjugate gradient algorithm, where it is 
necessary to store only the non zero entries of the matrix. 

The solution of equation (23), enables the computation of the electric field in the 
computational volume and the equivalent magentic and electric current densities on the 
surface terminating the computational domain. Using the equivalent electric and magnetic 
current densities on the surface terminating the computational domain, the radiated electric far 
field is computed as [4] 

17 / m ex P (-J k o r hc,zz — 

*W r) U~ = ~ jk ^o — Kr Jf(ee + <M>) 

y')expUk 0 sm(Q(x'costy + y sin(f>) + z'cosQ))dx'dy' 


+ A 


e x P( -jk a r) 
4nr 


r I* A A A A 

JJ (-6<j) + cf>0) 


M(*', y')exp(jk 0 sin cost}) + 'y sin<()) + z'cosQ))dx'dy' (24) 


where (r, 0, <(>) are the spherical coordinates of the observation point. The solution of 
equation (23) will also enables the calculation of electric field at the input plane, which can be 

used to calculte the reflection coefficient T at the input plane [9]. The input admittance is then 
calculated as 


( i-n 
(i + O ° 


(25) 


where Y 0 is the characteristic admittance of the feed transmission line. 
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Appendix 2 


Listing of the Distribution Disk 


/FEMOM3DR-1 . 0 


total 10 


drwxr-xr-x 

2 c j r 

1024 

Jul 

2 

09:47 

Examplel / 

drwxr-xr-x 

2 c j r 

512 

Jul 

2 

09:48 

Examp le2 / 

drwxr-xr-x 

2 c j r 

1024 

Jul 

2 

09:49 

Example3 / 

drwxr-xr-x 

2 c j r 

2048 

Jul 

2 

09:50 

FEM0M3DR/ 

drwxr-xr-x 

2 c j r 

512 

Jul 

2 

09:51 

PRE_FEM0M3 DR / 


/FEMOM3DR-1 . 0 / PRE_FEMOM3 DR 


total 63 

-rw-r— r— 

1 c j r 

202 Jul 2 

09:52 README 

-rw-r- -r- - 

1 cjr 

6741 Oct 31 

1997 cosmos2fem.f 

-rw-r— r— 

■ 1 cjr 

5295 Oct 31 

1997 edge . f 

-rw-r— r— 

1 cjr 

307 Oct 31 

1997 makefile 

-rw-r— r — 

1 cjr 

1076 Oct 31 

1997 meshin . f 

-rw-r— r-- 

1 cjr 

1723 Oct 31 

1997 paramO 

-rw-r— r— 

1 cjr 

1289 Oct 31 

1997 pmax . f 

-rw-r— r— 

1 cjr 

9 090 May 2 0 

14:02 pre_f emom3dr . f 

-rw-r— r — 

1 cjr 

3854 Oct 31 

1997 surf el . f 

/ F EM0M3 DR - 1 

. 0/FEMOM3DR 



total 380 

-rw-r — r— 

1 cjr 

472 Jul 

2 09.: 51 README 

-rw-r — r— 

1 cjr 

4825 Oct 

30 1997 analy . f 

-rw-r— r— 

1 cjr 

4195 Oct 

30 1997 basis. f 

-rw-r— r — 

1 c j r 

885 Oct 

30 1997 bess j . f 

-rw-r— r— 

1 cjr 

984 Oct 

3 0 1997 bess j 0 . f 

-rw-r— r— 

1 cjr 

1010 Oct 

30 1997 bessjl.f 

-rw-r— r — 

1 cjr 

3921 Jun 

11 09:40 bicgdns . f 

-rw-r— r— 

1 cjr 

2027 Oct 

30 1997 elembd.f 

-rw-r — r— 

1 cjr 

4090 May 

14 16:03 elmatr.f 

-rw-r--r-- 

1 cjr 

21317 Jul 

1 10:14 femom3dr.f 

-rw-r— r — 

1 cjr 

1887 Oct 

30 1997 f ourier_rwg . f 

-rw-r— r — 

1 cjr 

2640 Oct 

30 1997 fourierxy.f 

-rw-r— r — 

1 cjr 

658 Jun 

5 14:37 makefile 

-rw-r— r — 

1 cjr 

4405 Jul 

1 09:59 param 
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-rw-r--r-- 

-rw-r--r-- 

-rw-r—r — 

-rw-r--r-- 

-rw-r— r— 

-rw-r--r-- 

-rw-r— r— 

-rw-r--r-- 

-rw-r--r-- 

-rw-r--r-- 

-rw-r--r-- 

-rw-r--r— 

-rw~r--r— 

-rw-r—r— 

-rw-r— r— 

-rw-r--r— 

-rw-r--r — 

-rw-r--r-- 

-rw-r--r— 

-rw-r— r-- 

-rw-r--r-- 

-rw-r--r— 

-rw-r— r-- 

-rw-r— r-- 

-rw-r — r— 

-rw-r— r — 

/FEM0M3DR-1 

total 1684 
-rw-r--r— 
-rw-r--r— 
-rw-r— r— 
-rw-r— r-- 
-rw-r— r-- 
-rw-r— r-- 
-rw-r--r— 
-rw-r--r-- 
-rw-r--r — 
-rw-r— r — 
-rw-r--r-- 
-rw-r— r— 


1 

c j r 

882 

Oct 

30 

1997 

pleq . f 

1 

cjr 

4933 

Oct 

30 

1997 

quadpts . f 

1 

c j r 

3102 

May 

20 

14:07 

radpattn. f 

1 

cjr 

2983 Jun 2 

! 09:55 

> scatter_coax . 

1 

c j r 

4240 

May 

27 

16:29 

scatter_cwg . f 

1 

cjr 

4384 

Jul 

1 

10:32 

scatter_rwg. f 

1 

cjr 

307 

Oct 

30 

1997 

second. f 

1 

c j r 

4338 

Jun 

15 

09:45 

selmts__coax . f 

1 

c j r 

5312 

Oct 

30 

1997 

selmts_cwg . f 

1 

cjr 

9097 

Oct 

30 

1997 

selmts_rwg . f 

1 

cjr 

2219 

Oct 

30 

1997 

triangl_rwg. f 

1 

cjr 

2682 

Oct 

30 

1997 

triang_coax. f 

1 

c j r 

2611 

Oct 

30 

1997 

triang_cwg. f 

1 

c j r 

819 

Oct 

30 

1997 

triang_rwg . f 

1 

c j r 

1438 

Oct 

30 

1997 

triangeh . f 

1 

cjr 

2905 

Oct 

30 

1997 

triangej . f 

1 

cjr 

3089 

Oct 

30 

1997 

triange j 0 . f 

1 

c j r 

2449 

Oct 

30 

1997 

triangej 01 . f 

1 

cjr 

3105 

Oct 

30 

1997 

triangem. f 

1 

c j r 

1693 

Oct 

30 

1997 

triangemO . f 

1 

c j r 

1238 

May 

8 

10:08 

unorm. f 

1 

cjr 

469 

Oct 

30 

1997 

vcross . f 

1 

c j r 

382 

Oct 

30 

1997 

vdot . f 

1 

c j r 

5148 

May 

19 

08:37 

zmatrixeh . f 

1 

c j r 

9146 

Jul 

1 

10:12 

zmatrixe j . f 

1 

cjr 

7738 

Jul 

1 

10:12 

zmatrixem. f 


. O/Examplel 


1 

c j r 

22 

Jun 

2 

08 : 55 

coax. MAT 

1 

c j r 

75111 

Jun 

11 

09:31 

coax .MOD 

1 

cjr 

5050 

Jul 

1 

13:47 

coax . OUT 

1 

cjr 

133383 

Jul 

1 

13:15 

coax . PIN 

1 

c j r 

751 

Jul 

1 

13:15 

coax. POUT 

1 

c j r 

1113 

Jun 

11 

09:30 

coax . SES 

1 

cjr 

80 

Jul 

1 

13:47 

coax_bicgd . DAT 

1 

cjr 

317686 

Jul 

1 

13:15 

coax_edges . DAT 

1 

c j r 

93462 

Jul 

1 

13:15 

coax_nodal . DAT 

1 

c j r 

220044 

Jul 

1 

13:15 

coax_surfed.DAT 

1 

cjr 

12126 

Jul 

1 

13:15 

coax_sur f el . DAT 

1 

c j r 

37 

Jun 

15 

09:46 

input 
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/FEM0M3DR-1 . 0/Example2 


total 2302 



-rw-r--r— 

1 

cjr 

-rw-r--r-- 

1 

c j r 

-rw-r--r-- 

1 

cjr 

-rw-r--r-- 

1 

c j r 

-rw-r— r-- 

1 

c j r 

-rw-r— r — 

1 

cjr 

-rw-r— r— 

1 

cjr 

-rw-r— r— 

1 

c j r 

-rw-r— r-- 

1 

cjr 

-rw-r—r — 

1 

c j r 

-rw-r— r— 

1 

c j r 


/FEM0M3DR-1 . 0 / Examp le3 


32 

May 

21 

09:42 

input 

22 

May 

12 

15:24 

rwg . MAT 

106939 

May 

12 

15:23 

rwg . MOD 

4895 

Jul 

1 

10:30 

rwg . OUT 

187701 

Jul 

1 

09:50 

rwg . PIN 

751 

Jul 

1 

09 : 50 

rwg . POUT 

80 

Jul 

1 

10:30 

rwg_bicgd . DAT 

460142 

Jul 

1 

09:50 

^wg_edges . DAT 

134680 

Jul 

1 

09:50 

rwg_nodal . DAT 

263284 

Jul 

1 

09:50 

rwg__sur f ed . DAT 


16644 Jul 1 09:50 rwg_surfel.DAT 


total 

-rw-r- 

2475 

-r-- 

1 

c j r 

22 

-rw-r- 

-r-- 

1 

c j r 

115058 

-rw-r- 

-r-- 

1 

cjr 

4934 

-rw-r- 

-r-- 

1 

c j r 

202757 

-rw-r- 

-r- - 

1 

cjr 

751 

-rw-rw 

— 

1 

cjr 

1611 

-rw-r- 

-r-~ 

1 

cjr 

80 

-rw-r- 

-r-- 

1 

c j r 

506778 

-rw-r- 

-r-- 

1 

cjr 

147036 

-rw-r- 

-r-- 

1 

c j r 

265164 

-rw-r- 

-r-- 

1 

c j r 

19282 

-rw-r- 

-r~- 

1 

cjr 

32 


May 

21 

09:15 

cwg . MAT 


May 

27 

14:46 

cwg . MOD 


Jul 

1 

14:50 

cwg . OUT 


Jul 

1 

14 : 01 

cwg. PIN 


Jul 

1 

14:01 

cwg . POUT 


May 

27 

14:45 

cwg . SES 


Jul 

1 

14:50 

cwg_bicgd . 

DAT 

Jul 

1 

14 : 01 

cwg_edges . 

DAT 

Jul 

1 

14:01 

cwg_nodal . 

DAT 

Jul 

1 

14 : 01 

cwg_surf ed 

.DAT 

Jul 

1 

14:01 

cwg_surf el 

.DAT 

May 

27 

15:37 

input 
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Appendix 3 


Sample *.SES files of COSMOS/M 


The geometry modeling and meshing can be accomplished by using COSMOS/M. A 
variety of commands are available to define geometries. The constructed geometry is meshed 
and the mesh data can be written to a file with the Modinput command. Dielectric materials 
are identified by using material property command before meshing the corresponding part of 
the dielectric material. These are used as indices to tetrahedral elements, which will 
correspond to an entry in the problem. MAT file. Specification of the surfaces which are 
perfectly conducting, surfaces forming the radiating aperture and the input plane is 
accomplished by enforcing pressure boundary conditions on respective surfaces. Before the 
pressure condition is specified, a load condition has to be defined to indicate what type of 
surface is being specified. Load conditions of 1, 2, and 3 corresponds to perfectly conducting 
surface, surface at the fictitious outer boundary and surface at the input plane respectively. 

The *.SES files for the sample runs presented in section 4 are given below. 


Example 1: 

c* 

C* COSMOS/M Geostar VI. 75 

C;x Problem : /usrO /c j r/COSMOS/3d/FEMOM3DR/coax/coax 

Date : 

C* 

PLANE Z 0 1 
VIEW 0010 
PT 1 0 0 0 
PT 2 1 0 0 
SCALE 0 

CRPCIRCLE 1 1 2 1 360 4 

SCALE 0 

PT 6 1.57 0 0 

SCALE 0 

CRPCIRCLE 5 1 6 1.57 360 4 
SCALE 0 

CT 100. 5110 
CT 200. 5150 
RG 1 2 2 1 0 


35 



LT> 

CM 

CM 

I 

LD 

CM 

O 

LO 

CM 

CM 

LO 

CM 

CM 

LD 

CM 

O 

in 

CM 

CM 



CD 

no 


PRG 314 
CLS 1 
ACTSET LC 
PSF 521 



CLS 1 

ACTSET LC 3 
PRG 131134 

Example 2: 

c* 

C* COSMOS/M Geostar VI. 75 

C* Problem : /usrO /c j r/COSMOS/3d/FEMOM3DR/rwg 

7 

C* 

C* FILE rwg.in 1111 

SF4CORD 1 -0.35 -0.155 0 0.35 -0.155 0 0.35 0.155 0 -0 35 
0 

SCALE 0 

SFEXTR 1 4 1 Z -0.25 
SCALE 0 
CLS 1 

SF4CR 6 5 12 8 11 0 
PH 1 SF 1 0.1 0.001 1 
PT 9 -0.45 -0.25 0 
PT 10 0.45 -0.25 0 
PT 11 0.45 0.25 0 
PT 12 -0.45 0.25 0 
SCALE 0 

CRLINE 13 9 10 

CRLINE 14 10 11 

CRLINE 15 11 12 

CRLINE 16 12 9 

CT1 00.14123 40 

CT 200.14 13 14 15 16 0 

RG 1 2 2 1 0 

SFEXTR 13 16 1 Z -0.3 

CLS 1 

SF4CR 11 17 20 22 24 0 

SF4CORD 12 -0.5 -0.3 0.1 0.5 -0.3 0.1 0.5 0.3 0.1 -0503 

SCALE 0 ' 

SFEXTR 25 28 1 Z -0.5 
CLS 1 

SF4CR 17 29 36 32 35 0 
CLS 1 

PHPLOT 111 
SELINP SF 1 1 1 1 
SELINP SF 7 11 1 1 
SELINP RG 1 1 1 1 
CLS 1 


Date : 


0.155 


0.1 


37 


PH 2 SF 7 0.1 0.001 1 
CLS 1 

UNSELINP SF 1 1 1 1 
UNSELINP SF 7 11 1 1 
UNSELINP RG 1 1 1 1 
SELINP SF 12 17 1 1 
PH 3 SF 12 0.1 0.001 1 
PART 111 
PART 223 
CLS 1 

PARTPLOT 221 

PARTPLOT 121 

MPROP 1 PERMIT 1 

MA__PART 1110 0 4 

MA_PART 221004 

NMERGE 1 605 1 0.0001 000 

NCOMPRESS 1 605 

CLS 1 

INITSEL , SF , 1 , 1 
INITSEL , RG ,1,1 
CLS 1 

ACTSET LC 1 

PSF 2151114 

PSF 7 1 11 1 1 1 4 

PRG 111114 

ACTSET LC 2 

PSF 12 2 17 1 2 2 4 

ACTSET LC 3 

PSF 6361334 

Example 3: 

c* 

C* COSMOS/M Geostar VI. 75 

C* Problem : cwg Date : 5-27-98 Time : 15: 3:49 

C* 

C* FILE cwg. in 1111 
PLANE Z 0 1 
VIEW 0010 
PT 1 0 0 0 
PT 2 3.75 0 0 

CRPCIRCLE 1 1 2 3.75 360 4 
SCALE 0 

CT 101. 2110 
RG 1 1 1 0 
PT 6 -4.5 -4.5 0 
PT 7 4.5 -4.5 0 
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PT 8 4.5 4.5 0 
PT 9 -4.5 4.5 0 
SCALE 0 
CRLINE 567 
CRLINE 678 
CRLINE 789 
CRLINE 896 
CT 201. 2150 
RG 2 2 2 1 0 
SFEXTR 1 4 1 Z -3.0 
VIEW 1110 

RGGEN 1111000 -3.0 
PT 14 -4.5 -4.5 -3 .25 
PT 15 4.5 -4.5 -3.25 
PT 16 4.5 4.5 -3.25 
PT 17 -4.5 4.5 -3.25 
CLS 1 
CLS 1 
SCALE 0 
SCALE 0 

CRLINE 17 14 15 
CRLINE 18 15 16 
CRLINE 19 16 17 
CRLINE 20 17 14 
CT 401.21200 
RG 4 1 4 0 

SFEXTR 5 8 1 Z -3.25 
CLS 1 

SELINP SF 1 4 1 1 
SELINP RG 1 1 1 1 
SELINP RG 3 3 1 1 
CLS 1 

PH 1 SF 1 1.2 0.0001 1 
PART 111 
INITSEL SF 1 1 
INITSEL RG 1 1 
PT 18 -5 -5 0.5 
PT 19 5 -5 0.5 
SCALE 0 
PT 20 5 5 0.5 
PT 21 -5 5 0.5 
CRLINE 25 18 19 
CRLINE 26 19 20 
CRLINE 27 20 21 
CRLINE 28 21 18 
CT 501.21280 



RG 5 1 5 0 
RGGEN 15510 
SFEXTR 25 28 1 
CLS 1 

SELINP SF 5 12 
CLS 1 

SELINP RG 1 2 1 
SELINP RG 4 6 1 
CLS 1 

PH 2 RG 1 1.2 0 
PH 3 RG 5 1.2 0 
PART 223 
CLS 1 

INITSEL SF 1 1 
INITSEL RG 1 1 
PARTPLOT 121 
CLS 1 

PARTPLOT 111 
MA_PART 1110 
CLS 1 

PARTPLOT 221 
MA__PART 2 2 10 
PARTPLOT 121 
NMERGE 1 644 1 C 
NCOMPRESS 1 644 
ACTSET LC 1 
PSF 118111 
PRG 212114 
PRG 414114 
ACTSET LC 2 
PRG 526124 
PSF 9 2 12 1 2 2 
ACTSET LC 3 
PRG 333134 


0 0 -4.0 
Z -4.0 

1 1 

1 
1 

.0001 1 
.0001 1 

0 4 
0 4 

1.0001 0 0 0 
4 

4 
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Appendix 4 


Generic Input file format for PRE_FEMOM3DR 

The following is the format of the generic input file (problem.PIN) to be supplied to 
PRE_FEMOM3DR with required nodal data. 

• N n \ Number of nodes 

• N e : Number of trahedral elements 

• N p : Number of triangular elemets on PEC 
surfaces 

• N al : Number of triangular elements on 
surface at the outer boundary 

• N al : Number of triangular elements on 
surface at the input plane 

• N g : Maximum number of material groups 

x i>y»Zi 
x 2’ ^ 2 ’ z 2 

Coordinates of the nodes 


X N P ’ yN p ’ Z Np I 

n \V n 2\’ n 31’ n A\ m 8( 1) 

n 12’ n 22’ n 32’ n 42, mg(2) XT j 

Node numbers connecting each tetrahedral 

element 1, 2, 3, ,N e , and material group 

index number for each element 

n lN e ’ n 2N e ’ n 3N e ’ U 4N e , mg(N g ) 



41 



N e \ » n W n 2l’ n 3\ 

N e2 , n 12’ n 22’ n 32 

N eN p > n lN p ’ n 2N p> n 3N p 

^ e \ , n ll> n 2V n 3\ 
^ e 2, n \2’ n 22> n 32 


Global number of the terahedral element with a 
triangular face on PEC surface 

( N epN e2 , ,N eNp ) 

and three nodes connecting the triangular element 


Global number of the terahedral element with a 
triangular face on the outer boudary surface 

(^eV ^e2’ » N N ) 

a 1 

and three nodes connecting the triangular element 


N eN al ’ n lN al > n 2N al ’ n 3N ai 


N e \ * n W n 2V n 3\ 

N e2-> n l2’ n 22’ n 32 


Global number of the terahedral element with a 
triangular face on input plane surface 

(N e \, N e2 , , N e fsf a2 ) 

and three nodes connecting the triangular element 


N eN a2 , n lN a2 > n 2 N a2 ’ n 


3N n 
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